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The notion of fractality, in the context of positive-valued probability 
distributions, is conventionally associated with the class of Paretian prob- 
ability laws. In this research we show that the Paretian class is merely 
one out of six classes of probability laws - all equally entitled to be or- 
j^ ■ dained fractal, all possessing a characteristic power-law structure, and 

C"| all being the unique fixed points of renormalizations acting on the space 

of positive-valued probability distributions. These six fractal classes are 

""^J . further shown to be one-dimensional functional projections of underlying 

£^ ' fractal Poisson processes governed by: (i) a common elemental power-law 

O . structure; and, (ii) an intrinsic scale which can be either linear, harmonic, 

log-linear, or log-harmonic. This research provides a panoramic and com- 
prehensive view of fractal distributions, backed by a unified theory of their 
underlying Poissonian fractals. 
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1 Introduction 



Fractal objects are ubiquitous across many fields of science, and their study 
has attracted major interest by a broad array of researchers - see [I] -[I] and 
C^ ■ references therein. The geometric characteristic of fractals is invariance under 

changes of scale. The algebraic manifestation of scale-invariance is given by 
power-laws. 



'Department of Technology Management, Holon Institute of Technology, P.O. Box 305, 
Holon 58102, Israel. E-mail: eliazar@post.tau.ac.il 

t School of Chemistry, Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv 
69978, Israel. E-mail: klafter@post.tau.ac.il 



Power-laws facilitate the characterization of fractality when no geometry is 
present - the quintessential example being power-law probability distributions 
(see Chapter 38 in [J). 

Consider a population represented by a collection of points scattered ar- 
bitrarily on the positive half-line - the points representing the values of the 
population members. Examples include: earthquakes taking place in a given 
geological region, during a given period of time, measured by their magnitudes 
- each point representing the magnitude of an earthquake; stars in a given sector 
of space measured by their masses - each point representing the mass of a star; 
citizens of a given state measured by their wealth - each point representing the 
wealth of a citizen; insurance claims in a given insurance-portfolio measured by 
their costs - each point representing the cost of a claim; etc. 

Such populations are discrete objects possessing no natural geometry - and 
hence no natural geometric characterization of fractality. The natural setting 
for the analysis of such populations is statistical - providing the following con- 
ventional algebraic-statistic definition of fractality: a population is fractal if its 
population-values and their occurrence-frequencies are connected via a power- 
law. 

Shifting from the statistical perspective to the probabilistic perspective one 
picks at random a member of the population, and considers its random value X. 
Fractality, in the probabilistic setting, is characterized by a power-law survival 
probability of the random variable X: 

Vroh (x > x\ = (-Y (1) 

(x > a); the parameter a being an arbitrary positive lower bound, and the 
parameter a being an arbitrary positive exponent. 

The probability distribution corresponding to the survival probability of 
equation ([T]) is referred to as Paretian - named after the Italian economist 
Vilfrcdo Pareto who discovered, in 1896, a power-law distribution of wealth in 
human societies [5]. The Paretian probability distribution was empirically ob- 
served in a multitude of examples coming from diverse scientific fields [6]- [9] (see 
also the review [10] and references therein). 

The theoretical construction of the Paretian power-law probability distribu- 
tion is based on the following pair of foundations: (i) fractals are characterized, 
algebraically, by power-laws; (ii) probability distributions are characterized, 
statistically, by survival probabilities. 

The first foundation implicitly assumes that "fractality" is synonymous with 
"power-laws" . This implicit assumption is false. The notion of fractality - in 
the case of populations represented by arbitrarily-scattered real- valued points - 
can be defined from first principles. Namely, fractality can be defined via the 
elemental geometric notion of scale-invariance - rather than via the emergent 
algebraic notion of power-laws. This approach, undertaken in [11] . yields three 
classes of non-Paretian fractal populations. 



The second foundation implicitly assumes that random variables are uniquely 
characterized by their survival probabilities. This implicit assumption is, again, 
false. Indeed, there are many ways of characterizing a given probability distribu- 
tion (these characteristics will be rigorously defined in the sequel): Cumulative 
Distribution Functions (CDFs); Survival Distribution Functions (SDFs); Back- 
ward Hazard Rates (BHRs); Forward Hazard Rates (FHRs); Laplace Trans- 
forms (LTs); Moment Sequences (MSs); Log-Laplace Transforms (LLTs); Cu- 
mulant Sequences (CSs). 

Associating "fractal distributions" with power-law survival probabilities yields 
Paretian probability distributions. But what if we associate "fractal distribu- 
tions" with power-law Hazard Rates? or with power-law Log-Laplace trans- 
forms? or with power-law Cumulants? This question serves as the starting 
point of our research. 

This paper is devoted to the exploration of the definition of fractality in the 
context of positive-valued probability distributions. As we shall demonstrate, 
the notion of fractality is highly contingent on the distribution-characteristic 
used. Defining fractality via power-law structures of different distribution- 
characteristic leads to markedly different probability distributions including: 
Pareto, Beta, Frechet, Weibull, Levy Stable - all equally entitled to be consid- 
ered "fractal distributions" . 

Altogether we characterize six different classes of fractal distributions - 
each class emerging from a power-law structure of a different distribution- 
characteristic. Each of the six fractal classes characterized is shown to be associ- 
ated with a different renormalization: the members of each fractal class are the 
unique fixed points of a specific renormalization acting on the space of positive- 
valued probability distributions. Each of the six fractal classes characterized is 
also associated with a different Poissonian representation: the members of each 
fractal class are representable as a functional projection of an underlying class 
of Poisson processes defined on the positive half-line. 

Having characterized the six different classes of fractal distributions, their 
rcnormalizations, and their Poissonian representations, we turn to seek an un- 
derlying unifying fractal structure. To that end we study Poissonian renormal- 
izations - rcnormalizations of Poisson processes defined on the positive half- line 
- and characterize four classes of Poissonian fractals: the unique fixed points 
of multiplicative and power-law Poissonian rcnormalizations. 

The Poissonian fractals turn out to be governed by two structures: (i) a 
power-law structure common to all Poissonian fractal classes; (ii) an intrinsic 
scale which differentiates between the four Poissonian fractal classes and char- 
acterizes them. The intrinsic scale can be either linear, harmonic, log-linear, 
or log-harmonic. Moreover, the Poissonian fractals turn out to be the Pois- 
son processes underlying the aforementioned fractal distributions. And, the 
"algebraic fractality" on the "probability-distribution level" turns out to be a 
one-dimensional projection of a more elemental "geometric fractality" prevalent 
on the underlying "Poisson-process level" . 



This research provides a panoramic and comprehensive view of fractal dis- 
tributions, backed by a unified theory of their underlying Poissonian fractals. 
The manuscript is organized as follows. The six classes of fractal distributions 
- as well as their associated renormalizations - are characterized in Section O 
Poissonian representations of the fractal distributions are presented in Section 
[3J The Poissonian renormalizations and Poissonian fractals underlying the six 
classes of fractal distributions are unveiled in Section [4] 

Acronym glossary 

Throughout the manuscript the following acronyms shall be frequently used 
(the subsections in brackets indicate the location, in the manuscript, of the 
corresponding definitions): 

IID = Independent and Identically Distributed 

PDF = Probability Density Function (Subsection 12. ip 

CDF = Cumulative Distribution Function (Subsection 12. ip 

SDF = Survival Distribution Function (Subsection l2.ip 

BHR = Backward Hazard Rate (Subsection 12. 2p 

FHR = Forward Hazard Rate (Subsection 12. 2p 

LT = Laplace Transform (Subsection 12. 3p 

MS = Moment Sequence (Subsection 12. 3| 

LLT = Log-Laplace Transform (Subsection 12. 3| 

CS = Cumulant Sequence (Subsection 12. 3p 

CRF = Cumulative Rate Function (Subsection 13. ip 

SRF = Survival Rate Function (Subsection 13. ip 



The equality sign = shall henceforth denote equality in law (of random 
variables). 

2 Power-law characterization of fractal distribu- 
tions 

As noted in the introduction, positive- valued probability distributions have var- 
ious distribution-characteristics. In general, a distribution-characteristic Cd of 
a positive- valued probability distribution D is a function Cd — Cd{6) (S 6 9; 
being a subset of the non-negative half line) which uniquely determines D. 

The aim of this research is to explore the notion of fractality, in the context 
of positive-valued probability distributions, via the following definition: 



Definition 1 A probability distribution D is C-fractal if its distribution- characteristic 
Cd admits a power-law functional structure: 

C D (9)=c9~< (2) 

(9 G Q), where c is a positive coefficient and where 7 is a real exponent. 

Often, fractality is the manifestation of some underlying renormalization. In 
the context of positive-valued probability distributions a renormalization 1Z is 
a family of transformations 1Z = {TZ P } >0 - mapping probability distributions 
to probability distributions - which is consistent: A p- renormalization followed 
by a g-renormalization equals a pq-rcnormalization: 1Z P o lZ q — lZ pq (p, q > 0; 
the sign o denoting composition). 

A probability distribution D is a fixed point of the renormalization 1Z if it 
is a fixed point of each of the renormalization's transformations: 1Z P (D) = D 
(for all p > 0) . The connection between C-fractal probability distributions and 
renormalizations is given by the following definition: 

Definition 2 A renormalization 1Z is C-fractal if its set of fixed points coincides 
with the set of C-fractal probability distributions. 

In this section we study C-fractality with regard to each of the distribution- 
characteristics specified above. As shall be demonstrated, different distribution- 
characteristics will lead to very different meanings of fractality. 

2.1 Fractality via Frequencies 

Pareto's approach to analyzing the empirical data he gathered was based on 
frequencies: studying the occurrence- frequencies of the different population- 
values. In other words, Pareto focused on the Probability Density Function 
(PDF) /d(-) of a given probability distribution D. 

The PDF /d(-), in turn, induces the two most fundamental distribution- 
characteristics of a probability distribution D: (i) the Cumulative Distribution 
Function (CDF) F D (-), given by 

Fd(9) = I fo{x)dx (3) 

Jo 

(9 > 0); and, (ii) the Survival Distribution Function (SDF) Fd{-), given by 

/>oo 

F D {9) - / f{x)dx (4) 

Je 

(9 > 0). 

In this Subsection we study CDF-fractality and SDF-fractality. 



2.1.1 CDF-fractality 

The CDF Fjj(-) is monotone increasing from the level lim^o Fd{&) = to the 
level limg^oo Fjj(9) = 1. Hence, in order that the CDF Fpj(-) admit a power- 
law structure its underlying probability distribution D must be bounded from 
above. Admissible power-law CDFs are thus of the form 



F D (0) = {-) (5) 

(0 < 9 < a) , where the upper bound a and the exponent a are arbitrary positive 
parameters. 

With no loss of generality, the upper bound can be set to unity (a = 1) - 
yielding the Beta CDFs: 

F D (0)=9 a (6) 

(0<6»< 1). 

Let £ denote a random variable drawn from an arbitrary probability distri- 
bution D supported on the unit interval (0, 1). The conditional distribution of 
the scaled random variable £/p - contingent on the information that the random 
variable £ is no greater than the level p — is given by 



( K *M(^iS w 



(0 < p, 9 < 1). The conditional distribution of equation (J7J) induces the condi- 
tional renormalization 

F D { P 9) 

Fd( P ) 

{0<p,9< 1)Q 

A CDF Frj(-) is thus a renormalization fixed point if and only if it satisfies 
the functional equation Frj(xy) = Fd{x)Fd{v) (0 < x,y < 1). The solutions of 
this functional equation, in turn, are the Beta CDFs of equation ([6]). 

For probability distributions supported on the unit interval (0, 1) we conclude 
that: 

• The CDF-fractal probability distributions are the Beta distributions of 
equation ([6]). 

• The CDF-fractal renormalization is the conditional renormalization of 
equation ©. 



1 In this case the renormalization parameter p is restricted to the range < p < 1. 



2.1.2 SDF-fractality 

The SDF Fjj(-) is monotone decreasing from the level lim^o Frj(9) = 1 to the 
level limg^oo Fn(6) = 0. Hence, in order that the SDF Fr>{-) admit a power- 
law structure its underlying probability distribution D must be bounded from 
below. Admissible power-law SDFs are thus of the form 

F D {6) = Q a (9) 

(9 > a), where the lower bound a and the exponent a are arbitrary positive 
parameters. 

With no loss of generality the lower bound can be set to unity (a = 1) - 
yielding the Pareto SDFs: 

F D (9)=9- a (10) 

(0>1). 

Let £ denote a random variable drawn from an arbitrary probability dis- 
tribution D supported on the ray (l,oo). The conditional distribution of the 
scaled random variable £/p - contingent on the information that the random 
variable £ is greater than the level p — is given by 

Probfi> g |e>rU* rob «>^ (11) 

\p ' J Prob (£>p) ' 

(p,9 > 1). The conditional distribution of equation (fTTj) induces the conditional 
renormalization 

( Kp(Td))w = ^M) (12) 

(M>i)B_ 

A SDF Fd{') is thus a renormalization fixed point if and only if it satisfies 
the functional equation Fo(xy) = Fr)(x)F]j(jj) (x,y > 1). The solutions of 
this functional equation, in turn, are the Pareto SDFs of equation (|10j) . 

For probability distributions supported on the ray (1, oo) we conclude that: 

• The SDF-fractal probability distributions are the Pareto distributions of 
equation (JTUJ). 

• The SDF-fractal renormalization is the conditional renormalization of 
equation (fj"2j) . 



In this case the renormalization parameter p is restricted to the range p > 1. 



2.1.3 Exponential representations 

Both the aforementioned Beta and the Pareto probability distributions posses 
an underlying Exponential structure which we now describe. 

Let ^Beta denote a random variable governed by the Beta CDF of equa- 
tion ([5]); let £p arc to denote a random variable governed by the Pareto SDF of 
equation ([9]); and, let £ denote an Exponentially-distributed random variable 
with unit mean. It is straightforward to observe that the following exponential 
representations holds: 

£ B cta ^ a eX P { - - £ | aild Cparcto ^ « ex P \ ~ £ | ■ ( 13 ) 



Note that equation (113|) immediately implies a reciprocal connection between 
the Beta and the Pareto random variables: 

*. Law t Law J- f-t A\ 

SBcta — 7 and S Pareto — 7 ■ U 4 J 

^Pareto SBcta 



Let £ denote a random variable drawn from an arbitrary positive-valued 
probability distribution D. The conditional distribution of the translated ran- 
dom variable £ — p - contingent on the information that the random variable £ 
is greater than the level p - is given by 

p r o b (,- P >,u> P ) ^;° r b o «>;;;» 

(p, 8 > 0). The conditional distribution of equation (|15p induces the conditional 
renormalization 

J Mp- 

*d(p) 

(M>i). 

Equations (P~5|l - ifT^j) are the translational counterparts of equations 0-([5]| 
and equations pTjl -(fl2 |l . 

A SDF Fd(-) is a renormalization fixed point of equation (fT6| if and only 
if it satisfies the functional equation Fd(x + y) — Fd{x)Fd{v) {x,y > 0). 
The unique unit-mean solution of this functional equation is the unit-mean 
Exponential SDF. This characterizing property of the Exponential distribution 
- often referred to as "lack of memory" ([H]) Section XVII. 6) - is of prime 
importance in probability theory and its applications. As we see here, this 
elemental property also underlies the CDF-fractal renormalization and the SDF- 
fractal renormalization. 

2.2 Practality via Hazard Rates 

Let £ denote a random variable drawn from an arbitrary positive-valued prob- 
ability distribution D. What is the probability that the random variable £ be 



n p (F D ))(9) = ^ 7 ^i (16) 



realized at the level 9 - provided that it is not realized above the level 91 The 
answer to this question is given by the Backward Hazard Rate (BHR) Hd(-), 
defined as follows: 

(0>O). 

And what about the probability that the random variable £ be realized at 
the level 9 - provided that it is not realized below the level 91 The answer to 
this analogous question is given by the Forward Hazard Rate (FHR) H £>(■), 
defined as follows: 

*°W = £™o) p (t<e + 6\t>e) = i ^L (is) 

(9 > 0). The FHR plays a central role in Applied Probability and in the Theory 
of Reliability [13]- [11]. 

Both the BHR and the FHR are distribution-characteristics. Indeed, the 
CDF and the SDF can be reconstructed, respectively, from the BHR and the 
FHR via 

F D {9) =exp j- J H D (x)dx\ (19) 

(9 > 0), and via 

F D (9) = exp J - j H D {x)dx I (20) 

(0>O). 

In this subsection we study BHR-fractality and FHR-fractality. 

2.2.1 BHR-fractality 

As indicated above, the CDF Fu(-) is monotone increasing from the level 
lim^o Fd (9) = to the level lim^oo Fo(9) = 1. Hence, equation (fT9"]) im- 
plies that the BHR Hd(-) is integrable at infinity, and is non-integrable over 
the entire positive half-line (L Hi>(x)dx = oo). 

Admissible power-law BHRs thus yield the Frechet CDFs 

F D {9) = exp {-a9- a } (21) 

(9 > 0), where the coefficient a and the exponent a are arbitrary positive 
parameters. 

Let {£i, ■ ■ ■ ,£,„} denote a sequence of n IID random variables drawn from 
an arbitrary positive- valued probability distribution D. The distribution of the 
maximal random variable max {£ l5 • • ■ , £„} - scaled-down by the multiplicative 
factor n~ l ' a - is given by 

Prob (-1- max {£!,-• •,£«}<*) = (P™b (fi < n 1 /^))" (22) 



(9 > 0). The maximum distribution of equation (|22[) induces the maximal 
renormalization 



n p (F D )){9)=(F D ( P 1 l a 9)\ (23) 



(p,9>0). 

A CDF Fd(') is thus a renormalization fixed point if and only if its loga- 
rithm Gd{) — In (Fd(-)) satisfies the functional equation Go(xy) — x~ a Go{y) 
(x,y > 0). The solutions of this functional equation, in turn, are the Frechet 
CDFs of equation (j2"Tj) . 

For positive- valued probability distributions we conclude that: 

• The BHR-fractal probability distributions are the Frechet distributions of 
equation (|2"Tj) . 

• The BHR-fractal renormalization is the maximal renormalization of equa- 
tion ra. 



2.2.2 FHR-fractality 

As indicated above, the SDF Fd(-) is monotone decreasing from the level 
lime^o Fd{9) = 1 to the level \iu\e^oo Fr>(6) = 0. Hence, equation (J2T))) implies 
that the FHR Hjj(-) is integrable at the origin, and is non-integrable over the 
entire positive half- line (J Q Ho{x)dx = oo). 

Admissible power-law FHRs thus yield the Weibull SDFs 

F D (9)=cxp{-a9 a } (24) 

(9 > 0), where the coefficient a and the exponent a are arbitrary positive 
parameters. 

Let {^ , ■ ■ • , £ n } denote a sequence of n IID random variables drawn from 
an arbitrary positive- valued probability distribution D. The distribution of 
the minimal random variable minj^, • • ■ , £ n } — scaled-up by the multiplicative 
factor n 1 ' a - is given by 



Prob (n 1 /* ■ min {&, ■ ■ ■ , £„} > *) = (Prob U, > -^ j j (25) 

(9 > 0). The minimum distribution of equation (|25|) induces the minimal renor- 
malization 



(^(^))W=(^(^jJ ( 2 6) 



(M>o)._ 

A SDF Fd(-) is thus a renormalization fixed point if and only if its logarithm 
Gd(-) — In (Fd{-)) satisfies the functional equation Goi^y) = x a GD(y) {x, y > 
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0) . The solutions of the this functional equation, in turn, are the Weibull SDFs 
of equation (|24|) . 



For positive- valued probability distributions we conclude that: 

• The FHR-fractal probability distributions are the Weibull distributions of 
equation (|2"1)) . 

• The FHR-fractal renormalization is the minimal renormalization of equa- 
tion ((261) . 

2.2.3 Exponential representations 

Both the aforementioned Frechet and the Weibull probability distributions 
posses an underlying Exponential structure which we now describe. 

Let ^Frechet denote a random variable governed by the Frechet CDF of equa- 
tion (|2~Tjk let £ Weibull denote a random variable governed by the Weibull SDF of 
equation (|24[) ; and, let £ denote an Exponentially-distributed random variable 
with unit mean. It is straightforward to observe that the following exponential 
representations hold: 

Law (l \ 1/a Law (\ \~ 1/Q 

^Frechet = \~^J &n< ^ £ Weibull = Iq J ' (^7) 

Note that equation |27|) immediately implies a reciprocal connection between 
the Frechet and the Weibull random variables: 

£ Law t At L ™ ( r >a\ 

SFrechet — 7 anCL SWcibull — 7 ■ \ Z °> 

S Weibull S Frechet 

The Exponential distribution corresponds to the minimal renormalization of 
equation (|26[) with exponent a = 1. Indeed, the unique unit-mean solution of 
this renormalization is the unit-mean Exponential SDF. 

2.3 Practality via Laplace-space characteristics 

So forth, we considered fractality via "frequency-based" distribution-characteristics: 
CDFs, SDFs, BHRs, FHRs. In this Subsection we shift to Laplace space and 
turn to study fractality via the following "analytic-based" distribution-characteristics: 
Laplace Transforms; Moment Sequences; Log-Laplace Transforms; Cumulant 
Sequences. 

2.3.1 Laplace Transforms and Moment Sequences 

The Laplace Transform (LT) Ld(-) of a positive- valued probability distribution 
D is the Laplace Transform of its PDF /r>(') : 



/>CO 

L D {9)= / cxp{-to} 
Jo 



fo{x)dx (29) 
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(0 > 0). The LT is a distribution-characteristic - though, in general, the recon- 
struction of a PDF from a given LT is hard a task [16] . 

The LT Lu(-) is monotone decreasing from the level Lu(0) = 1 to the 
level lime^oo Lr>(9) = 0. Hence, LTs cannot admit the power-law structure of 
equation @. 

In case the LT Ld(-) admits a Taylor expansion around the origin, the 
Moment Sequence (MS) {Mrj(n)}^ =0 °^ the probability distribution D is well 
defined and is given by 

L D (9) = f]M D (n)t31 (30) 

n=Q 

(0 > 0). Reconstructing a probability distribution from a given MS is known as 
the Stieltjes Moment Problem [TO] . 

Since the Moment of order zero equals unity (Mu(0) = 1) MSs cannot admit 
the power-law structure of equation ([2]). 

We conclude that there are no LT-fractal and no MS-fractal probability 
distributions. 

2.3.2 Log-Laplace Transforms 

A "cousin" of the LT Ld(-) is its logarithm referred to as the Log-Laplace 
Transform (LLT) ^d(') and given by 

V D (9) = - hi (L D (9)) (31) 

(9 > 0). 

The LLT initiates at the origin (^d{0) = 0) and is monotone increasing 
(y' D {9) > 0) and concave (^>' D (9) < 0). 

Admissible power-law LLTs are thus the Levy Stable LLTs: 

# D (0)=a0 Q (32) 

ifi ^ 0); where a is an arbitrary positive coefficient and where the exponent a 
takes values in the range < a < 1. The Levy Stable LLTs of equation (|32|) 
admit the integral representation 



* D (e) = el~ex P{ -e x} (^- j ±;) dx 



(33) 



(0 > 0) - whose meaning will be explained in the sequel. 

(Apart from the special case a — 1/2, there is no "closed form" representa- 
tion for the PDFs of the Levy Stable probability distributions.) 
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Let {£ lt ■ ■ ■ , £„} denote a sequence of n IID random variables drawn from 
an arbitrary positive- valued probability distribution D. The LT of the aggregate 
£j + ■ • • + £ n - scaled-down by the multiplicative factor n -1 '" - is given by 

M- fk± ^})-(™{-^})" (34) 

(9 > 0). The LT of equation (|34f induces the aggregative renormalization 

(ll p (* D ))(6)=p* D (J^) (35) 

(p>0, 0> 0). 

A LLT &d{') is thus a renormalization fixed point if and only if it satisfies 
the functional equation ^ £,(xy) = x~ a ^D(y) (x,y > 0). The solutions of this 
functional equation, in turn, are the Levy Stable LLTs of equation (|32p . 

For positive- valued probability distributions we conclude that: 

• The LLT-fractal probability distributions are the Levy Stable distributions 
characterized by the LLTs of equation ([32]) . 

• The LLT-fractal renormalization is the aggregative renormalization of equa- 
tion (l35l). 



2.3.3 Cumulant Sequences 

In case the LLT ^d(-) admits a Taylor expansion around the origin, the Cumu- 
lant Sequence (CS) {Cd(«)}^ < L 1 of the probability distribution D is well defined 
and is given by 



*D(6) = -Y,C D (n)^- (36) 



(6 > 0). 

Power-law CSs of the form Coin) = an a (n = 1,2, •••), where the co- 
efficient a and the exponent a are arbitrary positive parameters, yield LLTs 
admitting the following integral representation: 



*d(o) = e£e xp {-ox} (f^Ta) (- ln ^))") dx 



(37) 



(0 > 0). 

The proof of equation |37|) is given in the Appendix; the meaning of this 
integral representation will be explained in the sequel. The renormalization 
associated with CS-fractal probability distributions is based on their underlying 
Poissonian structure - which, too, will be explained in the sequel. 

For positive- valued probability distributions we conclude that: 

13 



The CS-fractal probability distributions are characterized by the LLTs of 
equation (|37)) . 



• The CS-fractal renormalization is a Poissonian renormalization (yet to be 
presented) . 

2.4 Interim summery 

Table 1 summarizes the six classes of fractal probability distributions charac- 
terized in this Section. 

We note that in the context of IID sequences of positive- valued random vari- 
ables: (i) Extreme Value Theory asserts that the Frechet and Weibull distribu- 
tions are, respectively, the only possible linear scaling limits of the sequences' 
maxima and minima [T7]-[TH]; (ii) the Central Limit Theorem asserts that the 
one-sided Levy Stable distribution is the only possible linear scaling limit of the 
sequences' sums [2U]-[22]. 

Table 1 
The classes of fractal probability distributions 



Fractality 


Distribution 


Renormalization 


CDF-fractal 


Beta 


Conditional 


SDF-fractal 


Pareto 


Conditional 


BHR-fractal 


Frechet 


Maximal 


FHR-fractal 


Weibull 


Minimal 


LLT-fractal 


Levy Stable 


Aggregative 


CS-fractal 


— 


Poissonian 
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3 Poissonian representation of fractal distribu- 
tions 

In this section we provide Poissonian representations for all six classes of fractal 
probability distributions characterized in the previous Section. 

3.1 Poisson processes 

In this Subsection we recall the notion of Poisson processes. For further details 
the readers are referred to [25] . 

A Poisson process X on the positive half-line, with rate function r(-), is a 
random collection of positive-valued points satisfying the following properties: 
(i) the number of points Nx(I) residing in the interval / is a Poisson-distributed 
random variable with mean J I r(x)dx; and, (ii) if {Ik}^ is a finite collection of 
disjoint intervals then {Nx{Ik)\u is a finite collection of independent random 
variables. 

The rate function r(-) is the Poissonian analogue of the PDF in the context 
of probability distributions. The Poissonian analogues of the CDF and the SDF, 
respectively, are: (i) the Cumulative Rate Function (CRF) R(-), given by 



R(0) = / r {x)dx (38) 

Jo 

(9 > 0); and, (ii) the Survival Rate Function (SRF) R(-), given by 

/>oo 

R(9) = / r(x)dx (39) 



(0>O). 

The CRF R(-) is well defined if and only if the rate function r(-) is integrable 
at the origin - in which case it is a monotone non-decreasing function initiating 
from the origin (R(0) — 0). The SRF R(-) is well defined if and only if the 
rate function r(-) is integrable at infinity - in which case it is a monotone non- 
increasing function decreasing to zero (limg^oo R{9) = 0). 

The average number of points of the Poisson process X residing below the 
level is given by the CRF value R(9); the average number of points residing 
above the level 9 is given by the SRF value R(9). 

3.2 Poissonian maxima 

Consider the maximum of the Poisson process X, defined as follows: 

Xmax = max {x} . (40) 

xGPc 

We refer to the probability distribution of the random variable X max as the 
maximal distribution of the Poisson process X. 
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The maximum X max is smaller than the level 9 (9 > 0) if and only if the 
process X has no points residing above this level. Namely: {X max < 9} = 
{Nx((0, oo)) = 0}. Since the random variable Nx((0, oo)) is Poisson-distributed 
with mean R(9) , we obtain that the maximal distribution of the Poisson process 
X is characterized by the CDF 

F max (6)=exv{-R(9)} (41) 

(0>O). 

Equation (|41[) implies a one-to-one correspondence between Poisson pro- 
cesses (characterized by their SRFs R(-)) and their associated maximal distri- 
butions (characterized by their CDFs F max (-)). This one-to-one correspondence 
yields the following Poissonian representation of the CDF-fractal and BHR- 
fractal probability distributions: 

• A probability distribution D is CDF-fractal if and only if it is the maximal 
distribution of a Poisson process X with logarithmic SRF of the form 



R ^=- nw ) (42) 

(0 < 9 < a), where the upper bound a and the exponent a are arbitrary 
positive parameters. 

• A probability distribution D is BHR-fractal if and only if it is the maximal 
distribution of a Poisson process X with power-law SRF of the form 

R{9) = a9~ a (43) 

(9 > 0), where the coefficient a and the exponent a are arbitrary positive 
parameters. 

Equations ([4^]) and (|43[) follow, respectively, from equations (O and (f!H]) . 

3.3 Poissonian minima 

Consider the minimum of the Poisson process X, defined as follows: 

X min = min {x} . (44) 

We refer to the probability distribution of the random variable A m j n as the 
minimal distribution of the Poisson process X. 

The minimum X mm is larger than the level 9 (9 > 0) if and only if Poisson 
process X has no points residing below this level. Namely: {X mul > 9} — 
{Nx((0, 8]) — 0}. Since the random variable Nx((0,6]) is Poisson-distributed 
with mean R(9), we obtain that the minimal distribution of the Poisson process 
X is characterized by the CDF 

F min (9)=exp{-R(9)} (45) 
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(e > o). 

Equation (I45|) implies a one-to-one correspondence between Poisson pro- 
cesses (characterized by their CRFs R(-)) and their associated minimal distri- 
butions (characterized by their SDFs F m ; n (-)). This one-to-one correspondence 
yields the following Poissonian representation of the SDF-fractal and FHR- 
fractal probability distributions: 

• A probability distribution D is SDF-fractal if and only if it is the minimal 
distribution of a Poisson process X with logarithmic CRF of the form 

R(6) = -ln(^Y) (46) 

(9 > a), where the lower bound a and the exponent a are arbitrary positive 
parameters. 

• A probability distribution D is FHR-fractal if and only if it is the minimal 
distribution of a Poisson process X with power-law CRF of the form 

R(0) = a6 a (47) 

(6 > 0), where the coefficient a and the exponent a are arbitrary positive 
parameters. 

Equations (|4"6"|) and (|47[) follow, respectively, from equations © and ([M)l . 

3.4 Poissonian aggregates 

Consider the aggregate of the Poisson process X ', defined as follows: 

X agg = £ x . (48) 

We refer to the probability distribution of the random variable A agg as the 
aggregate distribution of the Poisson process X. 

The aggregate of equation (|4"8|) can be either convergent (A agg < oc) or 
divergent (A agg = oo). Campbell's theorem of the theory of Poisson processes 
(|23j. Section 3.2) implies that the aggregate is convergent if and only if the 
SRF R(-) is integrable at the origin - in which case the aggregate distribution 
of the Poisson process X is characterized by the LLT 

/•OO 

V am (0)=0 exp{-8x}R{x)dx (49) 

(0>o). 

Equation (|49[) implies a one-to-one correspondence between Poisson pro- 
cesses (characterized by their SRFs R(-)) and their associated aggregate distri- 
butions (characterized by their LLTs \t agg (-) - which, in turn, are characterized 
by the Laplace transforms of the underlying SRFs). This one-to-one corre- 
spondence yields the following Poissonian representation of the LLT-fractal and 
CS-fractal probability distributions: 
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• A probability distribution D is LLT-fractal if and only if it is the aggregate 
distribution of a Poisson process X with power-law SRF of the form 

^ = W^)W (50) 

(9 > 0), where a is an arbitrary positive coefficient and where the exponent 
a takes values in the range < a < 1. 

• A probability distribution D is CS-fractal if and only if it is the aggregate 
distribution of a Poisson process X with logarithmic SRF of the form 



R W = WTaj(- He) ) (51) 

(0 < 8 < 1), where the coefficient a and the exponent a are arbitrary 
positive parameters. 

Equations (|5D|) and (|5ip follow, respectively, from equations (|55|) and (|57|) . 

4 The underlying Poissonian fractals 

In Section [2] we characterized six classes of fractal probability distributions - 
each stemming from a different distribution-characteristic, and each associated 
with a different renormalization. In section [3] we have further seen that all six 
classes of fractal probability distributions admit Poissonian representations - 
either maximal, minimal, or aggregative. 

Is there any kind of an underlying order to this "little zoo " of fractal distri- 
butions? 

The answer is affirmative: all fractal distributions obtained are functional 
projections of underlying Poissonian fractals - as we shall show in this Section. 

4.1 Poissonian renormalizations and their fixed points 

In this Subsection we study renormalizations of Poisson processes defined on 
the positive half- line. We follow the renormalization approach used in |llj . 

4.1.1 Poissonian renormalizations 

Let \4> p } be a family of consistent scaling functions: monotone-increasing 
functions which map the positive half-line (0, oo) onto itself, and which sat- 
isfy the "consistency condition" (f> o <j) = <j) (p 1 q > 0; the sign o denoting 
composition). 

Given a Poisson process X with rate function r(-) we construct its p-order 
renormalization X p via the following two-step algorithm: (i) replace the process 
X by an intermediate Poisson process <Y' nt with rate function r' nt (-) = p ■ 
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r()|fj (ii) shift the points of the intermediate process X p nt using the p th scaling 
function <f> . The resulting p-order renormalization is given by 

x p = W*)}.^ • 

(The "consistency condition" is required in order to ensure that the Pois- 
sonian renormalization is consistent. Namely, that a p-order renormalization 
followed by a g-order renormalization equals a pg-order renormalization.) 

The connection between the CRF R p {-) and the SRF R p {-) of the p-order 
renormalization X p , and the CRF R(-) and the SRF R(-) of the original process 
X, is given by: R p {) — pR {4> p l {■)) andi? p (-) = pR (0 p 1 (-))i where the function 
</)~ 1 (-) denotes the inverse of the scaling function (f> p (-) (these results are an 
immediate consequence of the "displacement theorem" of the theory of Poisson 
processes - see Section 5.5 in [23]). 

Denoting by 1Z = {1Z P } >0 the Poissonian renormalization defined, we have: 

TZ p (R)=p(Ro ( f,- 1 ) and K p (~R) = p (R o 4> p x ) (52) 

(p>0). 

A Poisson process X is a fixed point of the renormalization 1Z if it is left 
statistically unchanged by the renormalization's action: the p-order renormal- 
ization X p being equal, in law, to the original process X. In terms of the CRF 
and SRF it is required that 1Z P (R) = R and 1Z P (i?) = R (for all p > 0). Using 
equation (|52|) we conclude that: the Poisson process X is a renormalization fixed 
point if and only if its CRF R(-) and SRF R{-) satisfy 

R o (j) p = pR and Ro <f) p = pR (53) 

(p>0). 

The two most fundamental Poissonian renormalizations are multiplicative 
and power-law. We now turn to characterize the fixed points of these renormal- 
izations. 

4.1.2 Fixed points of multiplicative renormalizations 

A multiplicative renormalization is based on a set of multiplicative scaling func- 
tions \4> p } ^ n - The consistency condition implies that the multiplicative scaling 
functions admit the form 

<f> p (x) = p 6 x (54) 

(x > 0), where the exponent e is an arbitrary non-zero parameter. 

The fixed points of a multiplicative renormalization governed by the scaling 
functions of equation ([54]) are as follows (the coefficient c being an arbitrary 
positive parameter): 



d Note that if p is integer then the intermediate Poisson process X* nt is the union of p IID 
copies of the original Poisson process X . 
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• In the case of a positive exponent e the renormalization fixed points are 
characterized by the CRF 

R{9) = c8 l/e (9 > 0). (55) 

The accumulation point of these fixed-point Poisson processes is a;* = oo|j 

• In the case of a negative exponent e the renormalization fixed points are 
characterized by the SRF 

R{9) = c9 1/e {9 > 0). (56) 

The accumulation point of these fixed-point Poisson processes is x* =0. 

4.1.3 Fixed points of power-law renormalizations 

A power-law renormalization is based on a set of power-law scaling functions 
\4> p ] ^ n - The consistency condition implies that the power-law scaling functions 
admit the form 

<j> p (x) = x pe (57) 

(x > 0), where the exponent e is an arbitrary non-zero parameter. 

The fixed points of a power-law renormalization governed by the scaling 
functions of equation (|57|) cannot range over the entire positive half-line (0, oo). 
Rather, they may range either on the unit interval (0,1) or on the ray (1, oo) 
(note that the power-law scaling functions of equation (|57p indeed map the unit 
interval (0, 1) and the ray (1, oo), respectively, onto themselves). 

The fixed points of a power-law renormalization governed by the scaling 
functions of equation (|57[) . and ranging over the unit interval (0,1), are as 
follows (the coefficient c being an arbitrary positive parameter) : 

• In the case of a positive exponent e the renormalization fixed points are 
characterized by the SRF 

R(8) = c{-\n9) 1/£ (O<0<1). (58) 

The accumulation point of these fixed-point Poisson processes is a;* =0. 

• In the case of a negative exponent e the renormalization fixed points are 
characterized by the CRF 

R(9)=c(-\n9) 1/e (O<0<1). (59) 

The accumulation point of these fixed-point Poisson processes is z* = 1. 



4 A point i, (0 < x* < oo) is said to be an accumulation point of the Poisson process X 
if, with probability one, there are infinitely many points of X within any given neighborhood 
of x*. 



20 



The fixed points of a power-law renormalization governed by the scaling 
functions of equation (|57[) . and ranging over the ray (l,oo), are as follows (the 
coefficient c being an arbitrary positive parameter): 

• In the case of a positive exponent e the renormalization fixed points are 
characterized by the CRF 

R{6) = c (In 9) 1/e (0>1). (60) 

The accumulation point of these fixed-point Poisson processes is x* = oo. 

• In the case of a negative exponent e the renormalization fixed points are 
characterized by the SRF 

R(9) = C (In 9) 1/e (0>1). (61) 

The accumulation point of these fixed-point Poisson processes is x* = 1. 

4.2 Poissonian fractals 

In the previous Subsection we obtained six classes of renormalization fixed- 
point Poisson processes. Excluding the fixed-point processes whose accumula- 
tion point is x* — 1 (an interior point of the positive half-line), and considering 
the fixed-point processes whose accumulation point is either the origin x* = 
or infinity x* = oo (the boundaries of the positive half-line), we now define four 
classes of Poissonian fractals. The fractal distributions of Section [2] shall turn 
out to be one-dimensional functional projections - either maximal, minimal, or 
aggregative - of these underlying Poissonian fractals. 

4.2.1 Linear Poissonian fractals 

The class of linear Poissonian fractals comprises of all Poisson processes gov- 
erned by CRFs admitting the power-law form 

R(9) = a6 a (9 > 0), (62) 

where the coefficient a and the exponent a are arbitrary positive parameters. 
The members of this class are fixed points of multiplicative Poissonian renor- 
malizations. For this class: 

• The maximal distribution is degenerate: the maximum X max equals infin- 
ity with probability one. 

• The minimal distribution is the FHR-fractal Weibull distribution, charac- 
terized by the SDF _ 

F mln (0)=exp{-a0 a } (63) 

(9 > 0). 

• The aggregate distribution is degenerate: the aggregate X agg is infinite 
with probability one. 
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4.2.2 Harmonic Poissonian fractals 

The class of harmonic Poissonian fractals comprises of all Poisson processes 
governed by SRFs admitting the power-law form 

R(9) = a6- a (6 > 0), (64) 

where the coefficient a and the exponent a are arbitrary positive parameters. 
The members of this class are fixed points of multiplicative Poissonian renor- 
malizations. For this class: 

• The maximal distribution is the BHR-fractal Frechet distribution, charac- 
terized by the CDF 

iW0)=exp{-a0- Q } (65) 

(0>O). 

• The minimal distribution is degenerate: the minimum X m i n equals zero 
with probability one. 

• If the exponent a is in the range < a < 1 then the aggregate distribution 
is the LLT- fractal Levy Stable distribution, characterized by the LLT 

*a gg (0) = r(l - a)aO a (66) 

(0>O). 

• If the exponent a is in the range a > 1 then the aggregate distribution is 
degenerate: the aggregate X agg is infinite with probability one. 

4.2.3 Log-linear Poissonian fractals 

The class of log-linear Poissonian fractals comprises of all Poisson processes 
governed by CRFs admitting the logarithmic form 

R{6) = a(ln{0)) a {6 > 1), (67) 

where the coefficient a and the exponent a are arbitrary positive parameters. 
The members of this class are fixed points of power-law Poissonian renormal- 
izations. For this class: 

• The maximal distribution is degenerate: the maximum X max equals infin- 
ity with probability one. 

• The minimal distribution is characterized by the SDF 

F nlin (9)=exp{~a(\n(0)) a } (68) 

(0>1). 

• The aggregate distribution is degenerate: the aggregate X am is infinite 
with probability one. 



The SDF of equation (|68|) reduces to the Pareto SDF when setting the 
exponent value a to unity. Thus, we refer to the distribution corresponding to 
this SDF as Hyper Pareto. 
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4.2.4 Log-harmonic Poissonian fractals 

The class of log-harmonic Poissonian fractals comprises of all Poisson processes 
governed by SRFs admitting the logarithmic form 

R{6) =a(- m(6>)) Q (0 < 9 < 1), (69) 

where the coefficient a and the exponent a are arbitrary positive parameters. 
The members of this class are fixed points of power-law Poissonian renormal- 
izations. For this class: 

• The maximal distribution is characterized by the CDF 

F m ^(9)=exp{-a(-\n(6)) a } (70) 

(0<6»< 1). 

• The minimal distribution is degenerate: the minimum X min equals zero 
with probability one. 

• The aggregate distribution is the CS-fractal distribution, characterized by 
the CS 

^ i s T(l +a)a ,„. 

C am (n) = — (71) 

(»=1,2,.-.). 

The CDF of equation ([70)) reduces to the Beta CDF when setting the expo- 
nent value a to unity. Thus, we refer to the distribution corresponding to this 
CDF as Hyper Beta. 

4.3 Structural properties of Poissonian fractals 

In this Subsection we describe the structural properties of the four classes of 
Poissonian fractals presented in the previous Subsection. 

4.3.1 Power-law structure and intrinsic scales 

The CRFs of the linear and log-linear Poissonian fractals admit the power-law 
structure 

R(-) = a(s(-)) a , (72) 

where the coefficient a and the exponent a are arbitrary positive parameters, 
and where the function S(-) is the intrinsic scale of the class under consideration: 

• Linear scale S(9) = 9 in the case of linear Poissonian fractals (9 > 0). 

• Log-linear scale S(9) — ln(#) in the case of log-linear Poissonian fractals 
(0>1). 
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Analogously, the SRFs of harmonic and log-harmonic Poissonian fractals 
admit the power-law structure 

R(-) = a(s(-)Y , (73) 



where the coefficient a and the exponent a are arbitrary positive parameters, 
and where the function S(-) is the intrinsic scale of the class under consideration: 

• Harmonic scale S(6) = 9~ in the case of harmonic Poissonian fractals 
(0>O). 

• Log-harmonic scale S(9) = In (0~ ) in the case of log-harmonic Poissonian 
fractals (0 < 6 < 1). 



All four classes of Poissonian fractals share the common power-law structure 
y = ax a . What distinguishes one fractal class from another is the intrinsic scale 
- on which the power-law structure is composed. 

4.3.2 Order statistics and Exponential representations 

All four classes of Poissonian fractals posses an underlying Exponential structure 
which we now describe. Let {£ n } n= i denote an IID sequence of Exponentially- 
distributed random variables with unit mean. 

The points of linear and log-linear Poissonian fractals can be listed in an 
increasing order £ x < £ 2 < £ 3 < •••. The order statistics {Cnl^Lii m turn, 
admit the following exponential representations: 

• Linear Poissonian fractals: 

/ f -L- -J- f \ 1 / a 

Law / <--l I" — i- ,, \ ,_ .. 

(n = 1, 2, • • • ). Equation (|74[) is the infinite-dimensional counterpart of the 
one-dimensional exponential representations of the Frechet distribution 
given in equation (127j) . 



Log-linear Poissonian fractals: 



t n L =e X J( £l + --- + £ A 1/a ) (75) 



(n = 1, 2, • • •). Equation (|75|) is the infinite-dimensional generalization of 
the one-dimensional exponential representations of the Pareto distribution 
given in equation {Tc 



Analogously, the points of harmonic and log-harmonic Poissonian fractals 
can be listed in a decreasing order ^ > £ 2 > £3 > ■ ■ • • The order statistics 
{£ n } n=1 , in turn, admit the following exponential representations: 
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T (P. A L P \ ~ y l a 

Law / t^l 



Harmonic Poissonian fractals: 

L L = ( ~ ' - ' ~ ) (76) 

(n = 1,2, • • •). Equation (JTSJ) is the infinite-dimensional counterpart of the 
one-dimensional exponential representations of the Weibull distribution 
given in equation (|27p . 



Log-harmonic Poissonian fractals: 



£i + --- + £,^ 1/q 



p n ^ exp _ ( — !_- ) l (77) 



(n = 1, 2, • • •). Equation (|77p is the infinite-dimensional generalization of 
the one-dimensional exponential representations of the Beta distribution 
given in equation (fT3|) . 

Equations (|74")) - (jTT)) follow from the "displacement theorem" of the theory of 
Poisson processes ([23], Section 5.5), combined with the fact that the increasing 
sequence {£i + • • • + £n} n= i forms a standard unit-rate Poisson process. 

The reciprocal connection between equations (j?4")) and (|7S|) is the infinite- 
dimensional counterpart of the one-dimensional reciprocal connection between 
the Frechet and Weibull distributions given by equation (|28|) : the reciprocal 
connection between equations (|75p and (|77p is the infinite-dimensional gener- 
alization of the one-dimensional reciprocal connection between the Beta and 
Pareto distributions given by equation (fT4")) . 

4.3.3 Transforming between fractal classes 

It is possible to transform from an "input" Poissonian fractal X belonging to one 
fractal class to an "output" Poissonian fractal y belonging to another fractal 
class via a simple point-to-point mapping x i— » y = ip(x) - which transforms 
the points x of the "input" Poissonian fractal to the points y of the "output" 
Poissonian fractal. 

The point-to-point mappings are given in Table 2, which should be read as 
follows: in order to transform from the fractal class of row % to the fractal class of 
column j one has to apply the point-to-point mapping y — ip(x) appearing in cell 
(i,j) of the Table. The construction of these point-to-point mappings follows 
straightforwardly from the "displacement theorem" of the theory of Poisson 
processes ([23J, Section 5.5). 

Each class of Poissonian fractals has one degenerate extremal and one non- 
degenerate extremal - see Table 3 below. The point-to-point mappings of Table 
2 transform the degenerate extremals amongst themselves, and transform the 
non-degenerate extremals amongst themselves. The point-to-point mappings of 
Table 2 do not, however, transform the aggregates of the different Poissonian 
fractal classes to each other. 
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Table 2 
Point-to-point mappings of the Poissonian fractal classes 





Lin. 


Har. 


Log-lin. 


Log-har. 


Lin. 


X 


x- 1 


exp(a;) 


cxp(— x) 


Har. 


x- 1 


X 


cxp(a; _1 ) 


exp(— x _1 ) 


Log-lin. 


ln(x) 


(M.*))- 1 


X 


x- 1 


Log-har. 


— In (a:) 


-anMr 1 


x- 1 


X 



Table 3 
Extremals of the Poissonian fractal classes 



Fractal 
Class 


Degenerate 
extremals 


Non-degenerate 
extremals 


Linear 


Max = oo 


Min ~ Weibull 


Harmonic 


Min = 


Max ~ Frechet 


Log-linear 


Max = oo 


Min ~ Hyper Pareto 


Log- harmonic 


Min = 


Max ~ Hyper Beta 
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5 Conclusions 

In this paper we examined the definition of fractality in the context of positive- 
valued probability distributions. We followed the conventional approach of asso- 
ciating the notion of fractality with power-law structures - considering various 
distribution characteristics, rather than the survival probability alone. 

We proved the existence of no less than six different classes of fractal prob- 
ability distributions - all admitting a characteristic power-law structure, and 
all being the unique fixed-points of rcnormalizations acting on positive-valued 
probability distributions. Each class manifested a markedly different meaning 
of fractality. 

All fractal classes were further shown to admit an underlying Poissonian 
structure - each fractal distribution being a one-dimensional functional projec- 
tion of an underlying Poisson process. The underlying Poisson processes, in 
turn, are fractal objects - being the unique fixed-points of Poissonian rcnormal- 
izations. 

The notion of fractality on the one-dimensional "probability-distribution 
level" emanated from the notion of fractality on the infinite-dimensional "Poisson- 
process level" : 

On the "probability-distribution level" fractality was defined algebraically 
via power-law structures, and the connection between the different classes of 
fractal distributions was unclear. 

On the elemental "Poisson-process level", however, fractality was defined 
via the geometric notion of population-renormalization, and a unified picture 
of fractality was obtained: it became vividly clear how all classes of fractal 
distributions emerge from the underlying Poissonian fractals, how they connect 
to each other, and how the underlying Poissonian fractals connect to each other. 

We have seen that on the elemental "Poisson-process level" , fractals do admit 
a universal power-law structure, yet this structure is intertwined with another 
key structure: the intrinsic scale which can be either linear, harmonic, log-linear, 
or log-harmonic. Whereas the power-law structure is common to all Poissonian 
fractals, it is the intrinsic scale which differentiates between the Poissonian 
fractal classes and characterizes them. 

This research provides a panoramic and comprehensive view of fractal dis- 
tributions, backed by a unified theory of their underlying Poissonian fractals. 
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6 Appendix: proof of equation ( 1371) 

We compute the LLT ^ d{@) (& > 0) corresponding to the power-law CS Cjj(ri) 
arT a (n = 1,2, •••). 

Note that 

j.a—1 



a f 00 i Q_1 

C D (n) = — =a exp {-nt} —-dt (78) 

n a J T(a) 



/o l(o.) 

(n = 1,2, • • • ). Hence, substituting equation (|75|) into equation ([55]) gives 



- E^=i ( a Jo°° ex P i~ nt } T{^ dt ) { -ij. 

JO ^Z^n=l n! y V r ( Q ) 

/ °°(l- exp {-fl exp {-«}}) (rfe^ 1 ) 



(79) 



rft 



(0>O). 

Now: 



*pW=^ 00 (l-exp{-9e X p{-t}}) L±-t a -Adt (80) 

(using the change of variables u = exp {— £}) 

(using integration by parts) 

= /%«p{-fe} ( /' « (-M^))"' 1 ^ dx (82) 

io \y a r(«) u y 

On the other hand (using the change of variables t = — m(u)) we have 
1 a (-ln(u)) a_1 /" ln(x) 



a .<*-!-■ a 



■ r dw = / ^7-T* di = 7^ r (- ln(a;)) a . 

I x r(a) u i T(a) r(l + a) 

(83) 
Hence, substituting equation (|8"3")l into equation ([5^1) we conclude that 



* D (9) = 0^ exp{-fe} L ( ^ a) (-ln(*))°) <fa (84) 

(e > o). 
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